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Abstract 

Biodiversity has reached a critical state. In this context, stakeholders need indicators that both provide a synthetic 
view of the state of biodiversity and can be used as communication tools. Using river fishes as model, we developed 
community indicators that aim at integrating various components of biodiversity including interactions between 
species and ultimately the processes influencing ecosystem functions. We developed indices at the species level 
based on (i) the concept of specialization directly linked to the niche theory and (ii) the concept of originality 
measuring the overall degree of differences between a species and all other species in the same clade. Five major 
types of originality indices, based on phylogeny, habitat-linked and diet-linked morphology, life history traits, and 
ecological niche were analyzed. In a second step, we tested the relationship between all biodiversity indices and land 
use as a proxy of human pressures. Fish communities showed no significant temporal trend for most of these 
indices, but both originality indices based on diet- and habitat- linked morphology showed a significant increase 
through time. From a spatial point of view, all indices clearly singled out Corsica Island as having higher average 
originality and specialization. Finally, we observed that the originality index based on niche traits might be used as an 
informative biodiversity indicator because we showed it is sensitive to different land use classes along a landscape 
artificialization gradient. Moreover, its response remained unchanged over two other land use classifications at the 
global scale and also at the regional scale. 
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Introduction 

In 2002, the 188 countries that are signatories to the 
Convention on Biological Diversity committed themselves to 
"achieve by 2010 a significant reduction of the current rate of 
biodiversity loss at the global, regional and national level" [1,2]. 
Even though this target was not achieved (the new target is 
2020), research in the field of biodiversity indicators has been 
growing during the last decade [3,4]. 



As biodiversity is a complex object and subject, a first step 
for improving conservation plans is to build indices, which are 
intended to synthesize and simplify data in quantitative terms. 
Indices vary depending on the biological level quantified, i.e. 
from genes to biomes. Such a variety of biodiversity levels 
respond to the numerous ways of examining biodiversity, as 
defined by the Convention of Biological Diversity [5]. As indices 
quantify an aspect of biodiversity, they can become useful 
indicators if they tell us about the impact of human pressures 
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on biodiversity. Facing global changes, species responses are 
not uniform [6]. Although a few species are not negatively 
affected by human activity and are flourishing, many are 
declining or will become extinct in the next century [7]. In this 
sense, the evaluation of biodiversity needs to move away from 
a reliance on species lists and case-by-case approaches to 
give a more global picture of what happens for most species in 
an ecosystem. 

Up until ten years ago, all river ecosystem indicators were 
assessed on their hydromorphological, chemical, and biological 
characteristics (e.g. IBGN [8] or EPT [9]). Because of their 
ability to integrate environmental variability at different spatial 
scales, fish assemblages have been studied and new 
indicators of river ecosystems have been developed (e.g. IPR 
[10], [11]). Although these indicators encompass the relative 
importance of geographic, ecoregional, and local factors, they 
were developed using the reference condition of pristine 
ecosystems without human impacts. As Baker and King (2013) 
point out, there is a crucial need for new aquatic indicators 
based on other criteria than biotic indices or summary metrics 
(e.g. taxon richness, ordination scores) especially in 
assessment and management [12]. Here, we develop a new 
approach, specifically dedicated to evaluate different functional 
approaches of fish biodiversity. It is the first study to synthesize 
comparisons between a large variety of fish traits: life history 
traits, morphological traits linked to habitat or diet, habitat 
niche, and an integrative approach based on abiotic habitat 
specialization. 

Community indices consider upper biological organization 
levels beyond the species level. They take into account the 
relationship between species inside the community, sometimes 
explicitly, such as in trophic networks [13,14], and sometimes 
implicitly, such as in niche or habitat specialization approaches 
[15]. Even if all community indices are species-based, they 
incorporate more complexity than species indicators because 
of these interactive approaches. They thus correspond more 
closely to a primary objective for indicators, producing a 
synthetic representation of biodiversity. Indeed, these 
indicators should help tackle the problem of maintaining the 
entire community integrity despite global changes by providing 
decision makers with more accurate information about human 
impacts on a global scale. In this respect, they are closer to the 
steady-state perspective, a frequently mentioned policy 
objective. It is easy to address the functional facet of 
biodiversity in this way and quite popular nowadays in the 
ecosystem "services" context [16]. However, basic summary 
metrics at the community level lose valuable information and 
non-linear declines should be undetected with aggregate 
responses [17]. In function of the study objectives it may be 
important, especially in conservation, to analyze the dataset 
species by species (or see TITAN, [17]) and it is always helpful 
to carefully interpret the community results. 

Finally, the criterion to create an indicator is to build good 
communication tools that are easy to understand and friendly 
to use, adapted to the context and scale of needs [18,19]. 
Indicators provide information to fuel dialogue between 
different scientific disciplines and stakeholders involved in 
biodiversity conservation. However, indicators also try to reach 



new targets identified as extremely important for the 
preservation of biodiversity. Indicators in this case attempt to 
open a dialogue and convince people not already involved in 
conservation including policy makers (local to international), 
judges, industrials, and farmers. 

In this paper we aim to develop indices to better understand 
functional patterns in space and time of river fish communities 
and to evaluate their potential as biodiversity indicators for 
environmental policy makers. 

First, we quantified the spatial and temporal changes in 
composition of French fish communities with two different 
approaches: originality and specialization indices. For the 
originality indices, we used four sets of functional traits (habitat 
niche, life history, diet-linked morphology, and habitat-linked 
morphology) and the phylogeny to obtain five matrices for the 
twenty-six common fish species considered. We first used the 
metric of originality defined as the rarity of species traits to 
obtain scores for each species [20]. Thus, the whole 
contribution of species to trait community depends on its 
originality. More precisely, as integrative community-traits 
indices, we computed the average value of the originality score 
depending on the density of species locally present. The 
second approach was based on niche theory and species 
specialization such as it has been done for birds [21,22]. To 
carefully interpret the community results we also explored 
spatio-temporal analyses at the species level. We identified 
regions of low originality or specialization communities at the 
national scale and explored the temporal changes through 
nineteen years. For the first time, we explored potential 
congruent or mismatched patterns between different functional 
traits approaches. 

Next, we evaluated the link between these community 
indices and human pressures via land use. We used land use 
as our proxy because threats to global freshwater biodiversity 
are mainly due to industrial and agricultural impacts [23]. We 
tested the sensitivity of each of the six indices to human 
pressures using habitat modification data sets, and used these 
results to select biodiversity indicators. Finally, we discussed 
the choice of indicators selected by communication criteria to 
give a clear message for stakeholders and, especially in our 
case, for environmental policy makers. 

Materials and Methods 

Fish database 

We worked with the database of the French National Agency 
for Water and Aquatic Environments (Onema), which contains 
records of standardized electrofishing protocols performed 
between 1990 and 2009 during low-flow periods (May- 
October). Electrofishing is considered the most effective non- 
destructive sampling procedure for describing fish assemblage 
structure [24]. Sampling protocols were defined depending on 
river width and depth. Streams were sampled by wading 
(mostly two-pass removal), while fractional sampling strategies 
were undertaken in larger rivers. Since the implementation of 
the EU Water Framework Directive's surveillance monitoring, 
protocols follow the recommendations of the European 
Committee for Standardization [25]. To compare inter-annual 
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Table 1. List of the studied species. 



Latin Name 


English common name 


French common name 


Abramis brama 


Common Bream 


Breme commune 


Alburnoides bipunctatus 


Bleak 


Spirlin 


Alburnus alburnus 


Bleak 


Ablette 


Ameiurus melas 


Black bullhead 


Poisson chat 


Anguilla Anguilla 


European Eel 


Anguille 


Barbatula barbatula 


Stone loach 


Loche franche 


Barbus barbus 


Common Barbel 


Barbeau fluviatile 


Carassius spp. 


Crucian carp 


Carassin 


Chondrostoma nasus 


Common nase 


Hotu 


Cottus gobio 


European Bullhead 


Chabot 


Cyprinus carpio 


Common carp 


Carpe commune 


Esox lucius 


Northern Pike 


B rochet 


Gasterosteus aculeatus 


Three spines stickleback 


Epinoche 


Gobio gobio 


Gudgeon 


Goujon 


Gymnocephalus cernuus 


Ruffe 


Gremille 


Lepomis gibbosus 


Pumpkinseed sunfish 


Perche soleil 


Leuciscus leuciscus 


Dace 


Vandoise 


Perca fluviatilis 


European perch 


Perche 


Pboxinus phoxinus 


Minnow 


Vairon 


Pungitius pungitius 


Ninespine stickleback 


Epinochette 


Rutilus rutilus 


Common roach 


Gardon 


Salmo salar 


Atlantic salmon 


Saumon atlantique 


Scardinius erytbrophthalmus 


Common rudd 


Rotengle 


Squalius cepbalus 


European chub 


Chevaine 


Telestes soufia 


Souffia 


Blageon 


Tinea tinea 


Tench 


Tanche 


doi: 10.1371/journal.pone.0080968.t001 



densities, however, only surveys performed with the same 
sampling protocol were selected in the whole data set. Fishes 
were identified to species level, counted and then released. We 
worked with the 26 species for which trait data were available 
(around 90% of the total abundance catch) (Table 1 ). 

We extracted collecting events from Onema's fish database 
using two different criteria: 

(i) All sites regardless of temporal coverage, which yielded 5 
403 sites with 1-20 years of sampling and a total of 13 076 
sampling occasions (Dataset 1). 

(ii) Only sites with at least 8 years of data, which yielded 557 
sites with 8-20 years of sampling from 1990 to 2009 and a total 
of 6942 sampling occasions (Dataset 2; see [26]). 

Trait dataset 

(i) : Habitat use. This dataset consists of five parameters 
describing the habitat use of a river species: foraging habitat, 
reproductive habitat, position in the water column, salinity 
tolerance and rheophily (the ability to live in fast moving water). 
The information has been gathered from different sources 
[27,28]. 

(ii) : Life history traits (LHT). The life history traits included 
in the study were: maximum lifespan, female age at maturity, 
number of spawns per year, logarithm of the maximum body 



length, relative fecundity (Number of ovocytes per gram), egg 
diameter, and parental care [27-29]. 

(iii) : Morphological database. Fourteen morphological 
traits related to two different axes of the niche (diet and habitat) 
were used (Figure 1, Table 2) [30-32]. Traits were measured 
from pictures collected mainly from FishBase [33,34], but for 
details see 35. All traits were standardized to account for 
differently sized photographs and species (e.g. standard 
length). 

(iv) : Phylogenetic dataset. We retrieved molecular data 
from three mitochondrial genes from GenBank (cytochrome b, 
cytochrome oxidase I and ribosomial 16S sub-unit). We 
inferred the best evolutionary model for each gene using 
maximum likelihood methods implemented in Paup4 ob10 [36]. 
The best model of molecular evolution was obtained using 
Modeltest based on the AIC criterion [37], (for more details see 
38). 

Human Pressure dataset 

The dataset of human pressures was provided by the 
European land-cover database CORINE, which classifies 
landscape units larger than 25 ha into one of 44 classes [39] on 
the basis of satellite digital images (e.g. SPOT and LANDSAT). 
We used the 2000 update and considered three alternative 
groupings of seven habitat classes: (i) The CORINE Land 
Cover (CLC) yields 5 habitat classes: "Forest", "Meadow", 
"Farming", "Urban", and a "Mix" (i.e. a Mix between agricultural 
and urban habitats) (ii) The EUROWATER (a special variant of 
CLC for freshwater common to the European scale), yields 6 
habitat classes with the addition of the "Intensive Urban" 
habitat class, and (iii) The ONEMA land use classification (a 
special variant of CLC and EUROWATER for freshwater 
common to the national scale) yields 7 habitat classes with the 
addition of the "Intensive Farming" class. Only the latter two 
classifications are used to test the reproducibility of our 
indicator. 

Here, we consider land use classification as a gradient of 
landscape artificialization under human pressures. Land use is 
a common proxy for human pressures in terrestrial 
communities [6,22]. The link between land use and human 
pressures in river has been reviewed at a global scale [40], but 
also on regional scale in North America [41,42] and Europe 
[43]. Marzin and al. (2013) showed a clear link between the 
CORINE Land Cover (CLC) dataset and different pollution and 
physical modifications at both local and regional scales [43]. If 
both human pressures are correlated with CLC, water quality 
parameters are more strongly correlated to land use than 
physical modifications. 

Statistical Analyses 

All statistical analyses were performed using R 2.15.1 (R 
Development Core Team. 2012), and more particularly the 
ade4 and nlme packages [44,45]. We calculated one index for 
each kind of data set, thus for a total of 6 indices: 4 functional 
originality indices, 1 phylogenetic originality index, and 1 
specialization index. 

(i): Functional and phylogenetic originality. To 
characterize the functional originality of each species, we used 
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Figure 1. Functional character measurements, a) All measurements are standardised by the standard length. Caudal peduncle 
length is also standardised by body depth; caudal peduncle depth was only standardised by body depth; b) pectoral fin position: 
pectoral fin dorsal side distance divided by body depth at pectoral fin; c) eye position: eye-ventral side distance divided by body 
depth at the eye; mouth position: mouth-ventral side distance at the position of the eye divided by body depth at the position of the 
eye. 

doi: 10.1371/joumal.pone.0080968.g001 



the mean of a set of functional traits from the different datasets 
described above. For each dataset a distance matrix was 
created using the Gower's dissimilarity index to allow the 
treatment of various statistical types of variables when 
calculating distances [46]. A hierarchical clustering (the 
unweighted pair-group clustering method using arithmetic 
averages: UPGMA) of the distance matrix produced a 
functional dendrogram comprising all the species. For each 
functional tree and the phylogenetic tree we used the 
procedure of Pavoine et al. [20] to estimate the biological 
originality of each species using the quadratic entropy of Rao 
[47]. Branch lengths and tree topology are jointly taken into 
account in the calculation of this index of originality. We 



computed both the Equal-split index [48] and the QE-based 
index [20]. The Equal-split index is more influenced by unique 
traits (trait states observed in a single species) than rare traits 
(trait states shared by a few species), whereas the reverse is 
true for the QE-based index. However, as both indices yielded 
similar results, we retained the equal-split index only, which is 
subsequently referred to as Species Originality Index (SOI). 
When it was possible we explored the sensitivity of our SOI to 
the addition of species in the data set [see the File S1]. 

(ii): Species Specialization Index. Ideally, specialization 
should be measured as the multi-dimensional breadth of a 
species' ecological niche. An integrative index of habitat 
specialization (Species Specialization Index, SSI) was 
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Table 2. Description of functional traits related to the habitat 
and diet niche axes [30-32]. From Schleuter et al. [35]. 





Niche axis 


Morphological Traits 


Code 


Functional Interpretation 




Pectoral fin length 


PL_SL 


Maneuvering speed, habitat 
velocity 


Habitat 


Vertical position of 
pectoral fin 


PFP 


Turning capacity 




Body Depth 


BD 


Maneuvering, hydrodynamics 
in the habitat 


Swimming 


Caudal peduncle 
Length / BD 


CpD 


Swimming ability 


ability 


Caudal peduncle 
Depth / BD 


CpL 


Swimming ability 




Caudal peduncle length 


CL 


Swimming speed 




Eyes Position 


EP 


Vertical position in the Water 
column 








Adaptation (i) light (turbidity 




Eyes diameter 


EDSL 


and diurnal) (ii) Relative prey 
size 


Diet 


Mouth Position 


MP 


Location of food acquisition 


and 


Length of longest 
barbell 


BarL^SL 


Non visual food detection, 
benthic feeders 


Food 


Head length 


HLSL 


Relative prey size 


acquisition 


Length of upper jaw 


MS 


Relative prey size 








Actual prey size (in 




Maximum size 


Lmax 


combination with head and 
upper jaw length) 



From Schleuter et al. [35]. 

doi: 10.1371/journal.pone.0080968.t002 



developed for birds [21], as the coefficient of variation 
(standard deviation/average) in average density of a species 
across habitats. We tested the relevance of this index in fishes. 
Because ecological habitat classes were missing for several 
species, we used habitat traits and four abiotic variables: 
temperature (sum of January to June air temperatures), 
longitudinal gradient, log of elevation, and slope (see 13 for 
more details). We had to take into account the geographical 
bias in the data set. This bias was linked to an over-sampling of 
headwaters. We therefore reassigned all the sampled points 
into 7 habitat classes with an approximately equal amount of 
samples in each habitat class. 

(iii): Community Indices. Each species can be ranked 
along a continuous gradient from the least to the most original 
or specialized species (X,,..., X;). Any species assemblage at 
time f can be characterized by the average specialization or 
originality taken across all individuals in the assemblage. 
These community level indices are simple weighted averages, 
i.e. T( a i,tXiVT a <,t> where a u is the relative abundance of species / 
in the assemblage at time f and X, the originality or 
specialization of species /'. In the following, CSI, = £(a,,SS/,)/ 
£a ; , is the Community Specialization Index and COI = £ 
(a it SOiyZa it , the Community Originality Index at time t. We 
explored the temporal and spatial variation of both community 
indices, COI and CSI, using mixed-effects linear models with 
sampling site as a random effect [49,50] and Akaike 



Information Criteria (AIC) model selection. We also explored 
the link of all the community indices between them by exploring 
their correlations by performing linear model. For the statistical 
independence of the data spatio-temporal effects and their 
interactions were taken into account, and we selected the 
model in function of its AIC. 

(iv): Community Indicators. We tested the relationships 
between CSI and five COI (four functional and one 
phylogenetic originality indicators) and landscape variables 
using mixed-effects linear models with sampling site as a 
random effect. Temporal (year) and spatial effect (geographical 
coordinates and watersheds) with their interactions were also 
taken into account. Because no R-squared can be calculated 
with random effect, we only obtained a proxy of the R-squared 
with the same model without the random effect. We used the 
CORINE Land Cover dataset (Figure 2) and its two variations 
to evaluate the reproducibility of our results and thus the 
sensitivity of each community index through habitat 
classifications. Then we studied the scale dependence of the 
community index response by exploring the relationship at the 
regional watershed scale. 

Results 

(i). Species Originality and Specialization Indices 

The four trait distance matrices can be visualized using trees 
(Figure 3). Trees based on life-history traits (Figure 3a), 
functional niche (Figure 3b), and diet-linked morphological 
traits (Figure 3c) were well balanced in the sense of Blum et al. 
[51]. These authors defined the balance of a tree as the 
average balance of its nodes, "assuming that a given node is 
completely balanced if it splits the sample into two subsamples 
of equal size". At the opposite, the tree based on habitat-linked 
morphological traits (Figure 3d) was highly unbalanced by the 
European eel (Anguilla anguilla), and, to a lesser extent by the 
groups common bream (Abramis brama), crucian carp 
(Carassius sp.), ninespine stickleback (Pungitus pungitius), 
three spines stickleback (Gasterosteus aculeatus). Using the 
equal-split metric, we computed four originality indices for each 
species to evaluate the three functional datasets and the 
phylogeny (Figure 4). As expected, A. anguilla was 
characterized by a high originality score for the habitat-linked 
morphological trait matrix (SOI = 0.81 compared to a mean of 
0.24). The two other imbalanced nodes had smaller originalities 
(Abramis group = 0.46 and Pungitius group = 0. 37). The 
species specialization index ranked the European eel as the 
most specialist species and the common bleak, Alburnus 
alburnus as the most generalist species (Figure 4). At the 
community level, the habitat-linked morpho-COl was sensitive 
to the presence of the European eel and to a lesser extent to 
the presence of the common bream, Abramis brama, and the 
crucian carp, Carassius sp. We tested the sensitivity of the 
Species Originality Index (SOI) based on traits to the addition 
of species in the initial dataset. The life history traits index and 
the habitat-linked morphological index were strongly correlated 
(respectively R 2 =81, R 2 =86). The diet-linked morphological 
index and the niche index were less correlated (respectively R 2 
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200 km 




J Urban (artificial lands) j |] Meadow 
|U Farming (except meadow) | | Forest (semi-natural lands) 



Watershed 



Figure 2. Spatial distribution of human pressures. We used the CORINE Land Cover data set and its land use classification as 
an artificialization gradient from natural habitat to urban habitat. Colors from green to red represent increasing pressures. 

doi: 10.1371/journal.pone.0080968.g002 



= 65 and R 2 =68) and thus more sensitive to the addition offish 
species in the initial dataset [see the File S1]. 

(ii). Community Indices: spatial and temporal patterns 

All statistical models retained by the AIC, with both datasets, 
contained the same variables: geographic coordinates, year, 
and their interaction, except for the Diet-COl where the 
watershed effect gave a better model (AIC). Because Corsica 
appeared to be an outlier (Figure 5), we re-ran all analyses 
excluding data from this area. With Corsica excluded, we found 
that watershed was a better spatial effect than geographic 
coordinates (AIC). 

Corsica clearly comes out as a hotspot of fish originality and 
specialization (Figure 5). In contrast, the Seine watershed 
presented the lowest originality and the most generalist fish 
communities. The CSI, habitat-linked morpho-COl and LHT- 
COI presented limited variation among sites and appeared to 
be ineffective to discriminate sites. In contrast, the diet-linked 



morpho -COI highlighted a strong originality in all small rivers, 
especially mountainous streams. Although the temporal effect 
was always retained in statistical models on the AIC, it was not 
always significant. However, Niche and both Morpho-COl 
increased significantly over the last years (Table 3). 

(iv). Community Indicators 

Our results showed that all COI, which are based on the 
same originality metrics, tended to be correlated with each 
other, while the CSI, based on a different approach and metric, 
was correlated with a subset of the COI only. All these 
correlations took into account space and time effect for the 
statistical independence of data. The habitat-linked morpho - 
COI was correlated with the LHT - COI (R 2 = 0.47, F 8 12829 = 
1407, P < 0.001), and with the diet-linked morpho -COI (R 2 = 
0.24, F 8 12829 = 5 1 8, P < 0.001). Interestingly, the phylogeny 
approach seemed to capture different proportion of 
morphological variation in function of the diet or the habitat 
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(a) 



1 r£ 




Salmo salar 
Anguilla anguilla 
Telestes soufia 
Leuciscus leuciscus 
Cottus gobio 
Ameiurus melas 
Perca fluviatilis 
Esox lucius 
Barbus barbus 
Abramis brama 
Scardinius erythrophtalmus 
Rutilus rutilus 
Chondrostoma nasus 
Pungitius pungitius 
Cyprinus carpio 
Squalius cephalus 
Carassius spp 
Tinea tinea 
Lepomis gibbosus 
Gasterosteus aculeatus 
Alburnoides bipunctatus 
Phoxinus phoxinus 
Alburnus alburnus 
Gymnocephalus cernuus 
Gobio gobio 
Barbatula barbatula 



(b) 




Gasterosteus aculeatus 
Pungitius pungitius 
Rutilus rutilus 

Scardinius erythrophtalmus 
Salmo salar 

Alburnoides bipunctatus 
Leuciscus leuciscus 
Alburnus albumus 
Squalius cephalus 
Lepomis gibbosus 
Telestes soufia 
Perca fluviatilis 
Esox lucius 
Phoxinus phoxinus 
Abramis brama 
Cyprinus carpio 
Gobio gobio 
Barbus barbus 
Chondrostoma nasus 
Anguilla anguilla 
Cottus gobio 
Ameiurus melas 
Carassius spp 
Tinea tinea 
Barbatula barbatula 
Gymnocephalus cernuus 



(c) 



(d) 





Gasterosteus aculeatus 
Pungitius pungitius 
Phoxinus phoxinus 
Gobio gobio 
Barbatula barbatula 
Cottus gobio 
Alburnoides bipunctatus 
Alburnus alburnus 
Telestes soufia 
Gymnocephalus cernuus 
Ameiurus melas 
Chondrostoma nasus 
Leuciscus leuciscus 
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Figure 3. Functional trees, (a) Tree based on Life History Traits (b) Tree based on functional Niche traits (c) Tree based on 
functional Diet Morphological traits (d) Tree based on functional habitat-linked morphological traits. 

doi: 10.1371/journal.pone.0080968.g003 



niche axis [phylo -COI and habitat-linked morpho -COI (R 2 = 
0.53, F 812829 = 1842, P < 0.001), and phylo-COl and diet-linked 
morpho COI (R 2 = 0.39, F 812829 = 1037, P < 0.001)]. Unlike other 
vertebrates such as birds or mammals (birds: [52]), the fish Life 
History Traits index presented a weak correlation with the 
phylo-COl (R 2 = 0.13, F 812829 = 239, P < 0.001). The niche-COl 
was more strongly correlated with the phylogenetic index (R 2 = 
0.56, F 812829 = 2013 10 4 , P < 0.001) than with the morphological 
indices (Diet: R 2 = 0.33, F 8 1282g = 798, P < 0.001; and Habitat: 
R 2 = 0.50, F 812829 = 1596, P < 0.001) even though the latter are 
assumed to represent niche axes. The niche-COl was also 



correlated with the LHT -COI (R 2 = 0.23, F 812829 = 475, P < 
0.001). 

The CSI was strongly correlated with the LHT-COI (R 2 = 0.74, 
F 8 ,12829 = 3-43 10 4 , P < 0.001) and habitat-linked morpho - COI 
(R 2 = 0.47, F 812829 1397, P < 0.001) but weakly with the niche 
(R 2 = 0.11, F 8 ,i2829 208, P < 0.001), and not with the diet-linked 
morpho (R 2 = 0.04, F 812829 65, P < 0.001) and phylo - COI (R 2 = 
0.08, F 812829 1 49, P < 0.001). It is important to note that the 
level of specialization measured here is more relevant to 
describe the Fish Life History Traits component than the habitat 
niche component. 
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Figure 4. Results of species values of originality and specialization in a Cleveland's dot plots in proportions (i.e. the 
sum of species values is one for each index). 

doi: 10.1371/joumal.pone.0080968.g004 



The sensitivity to human pressures of the six community 
indices was evaluated with respect to land use data and two 
variations of CORINE Land Cover (Table 4). All indices 
correlated with land use (Table 4-6), but with some variation. 
For example, some indices were sensitive to the different 
human pressures (farming or urban) represented here by an 
artificialization gradient (Figure 6). In contrast, the CSI was 
significantly higher for urban area than for agricultural or 
natural habitats (Figure 6). We used two variations of CORINE 
Land Cover (ONEMA and EUROWATER) to get an estimate of 
the community indices reproducibility in function of the arbitrary 
habitat classifications [Table 4-6, see the File S2] and only one 
COI was robust to the effect of habitat classifications: the Niche 
- COI (Table 4-6). The response of this latter index was also 
significantly sensitive to the different type of human pressures 
with a consistent behavior at national and regional scales 
(Figure 7). Within each watershed or over all watersheds the 
relationship between human pressures and Niche-COl is 
negative when it is expected to be negative (e.g. under human 
pressures like farmland and urban habitat) and reciprocally 
(e.g. under natural habitat). 



Discussion 

The aim of this study was to develop new functional indices 
for river fish and to evaluate their potential application as 
functional biodiversity indicators. For the first time in fish 
communities, we examine spatio-temporal patterns of six 
functional facets of biodiversity relying on two different 
theoretical approaches: specialization and originality. We 
identified common conservation priorities but also spatial 
mismatching in function of the trait considered. Then, we linked 
them with human land use pressures and we identified the 
community functional originality index based on niche traits as 
the most likely to become a functional biodiversity indicator. Its 
sensitivity to the nature and intensity of human disturbance, 
considered here by an artificialization gradient, at both regional 
and national scales, results in a simple message to 
communicate with policy makers and biodiversity managers. 

I. Community Indices 

There is a growing consensus that functional diversity based 
on species traits is a better predictor of ecosystem functioning 
than species number per se [53]. Species richness is currently 



PLOS ONE | www.plosone.org 



8 



November 2013 | Volume 8 | Issue 1 1 | e80968 



Freshwater Biodiversity Indicator 




Diet-COl _ Habitat - COI Phylogenetic - COI 




0.26-0.50 

0.51 -0.75 200 km 

• >0.76 



Figure 5. Spatial distribution of community specialization (CSI) and originality (COI) indices for 2005 (the years with the 
most sample stations). Colors from green to red represent increasing index values. 

doi: 10.1371/joumal.pone.0080968.g005 



the most used biodiversity index (and indicator) but it is highly 
scale-dependent, with local increases that are often 
accompanied by regional or global decreases and increases in 
between-site similarity [54]. Particularly, fish species richness 
tends to increase from upstream to downstream [55]; and the 
upstream part of many French rivers sustain only a few fish 
species (< 5 fish species, [56]). Species richness is thus an 
inadequate surrogate in the context of ecosystem function 
unlike community traits approaches, which appear more and 
more relevant in the literature especially to examine ecosystem 
integrity [16,22,57]. Community-trait indices take into account 
the species present in the area considered, species being 
grouped depending on their ecological or phylogenetic affinity. 
These indices compute the average value for a trait or 
character, depending on the frequency of species locally 
present. So, we did not consider intra-specific variability, which 
sometimes represents a significant proportion of the variability 
and complex spatio-temporal dynamic [58,59]. Even though 
specialization community indices seem to give the same 
message with presence/absence data [60], it is not always the 
case when a process or a function is measured using 



functional traits [19]. Moreover one of the most interesting 
points to use common species is based on the assumption that 
abundance plays an important role in ecosystem functioning 
[57,61]. 

Community Specialization Index (CSI) is a different approach 
than Community Originality Index (COI) approaches. If CSI is 
not clear on the underlying mechanisms explaining the precise 
ecosystem function, its well-known power comes from its 
holistic habitat approach. The central focus of CSI is not the 
species feature but its interaction with the environment by the 
habitat approach. This statement is closer to the Grinnell niche 
theory approach than the Hutchinson one [62]. And thus, in the 
CSI approach, the crucial point is the relevance of habitat 
description, not the species traits data set. 

On the other side, with the COI we study the distinctiveness 
of precise species traits and lineages, and thus we postulated 
that trait variation among species variation relates to functional 
differences in the ecosystem, which allows an interpretation in 
terms of ecosystem function or "services". The set of traits 
selected is a crucial step toward this goal, especially if we want 
the dynamics of the indices to reflect ecosystem function [57]. 
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Table 3. Results on temporal effect for three functional 
community indices, all the other indices present no 
temporal effect. 



Temporal main effect 




Datasetl 


Dataset2 


COI- Niche 


Coef. 


0.002 


0.002 


Year*(x + y) 


t value 


2.65 


2.32 




p value 


0.008 


0.02 


COI- Mhab 


Coef. 


0.004 


0.005 


Year*(x + y) 


t value 


3.67 


3.11 




p value 


O.001 


0.002 


COI- Mdiet 


Coef. 


0.005 


0.003 


Year*(x + y + w) 


t value 


3.45 


2,09 




p value 


O.001 


0.04 



The main temporal effect is given for each best model selected on the AIC, each 
statistical model has been run with a station as a random effect and as a fixed 
effect time and space and their interactions. The main difference is from the final 
space effect, which could be abscissa (x), ordinate (y), watershed (w) effect and all 
possible combinations. Results are similar without Corsica Island. 



doi: 10.1371/journal.pone.0080968.t003 

Table 4. Sensitivity of community indices to human 
pressures as defined and measured in CORINE Land 
Cover. 



Index 


Value 


Meadow 


Farming 


Mix 


Urban 


CSI 


Coef 


-0.04 


NS 


NS 


0.1 


R 2 = 0.26 


t-value 


-7 


NS 


NS 


8 




p-value 


<0.001 


NS 


NS 


<0.001 


LHT-COI 


Coef 


-0.03 


-0.01 


NS 


0.03 


R 2 = 0.28 


t-value 


-10 


-7 


NS 


4 




p-value 


<0.001 


O.001 


NS 


<0.001 


Niche-COl 


Coef 


-0.01 


-0.02 


NS 


-0.02 


R 2 = 0.27 


t-value 


-4 


-10 


NS 


-3 




p-value 


<0.001 


O.001 


NS 


0.04 


Diet-COl 


Coef 


-0.006 


-0.008 


NS 


-0.03 


R 2 = 0.06 


t-value 


-3 


-6 


NS 


-6 




p-value 


0.006 


<0.001 


NS 


<0.001 


Habitat-COl 


Coef 


NS 


NS 


NS 


0.06 


R 2 = 0.38 


t-value 


NS 


NS 


NS 


4 




p-value 


NS 


NS 


NS 


<0.001 


Phylogenetic-COl 


Coef 


NS 


-0.03 


NS 


-0.1 


R 2 = 0.21 


t-value 


NS 


-6 


NS 


-5 




p-value 


NS 


<0.001 


NS 


<0.001 



The relation between the community indices and the land used modification has 
been performed with a mixed-effects linear model with sampling site as a random 
effect and temporal (year) and spatial effect (geographical coordinates or 
watershed) and their interactions. The land use effects are given for each best 
model selected on the AIC. The coefficient effect of each habitat class is relative to 
the natural class (thus forest coefficient is always 0). We corrected the p-values for 
multiple tests using the Benjamini Yekutieli False Discovery Rate. We added the 
Rsquared adjusted from the equivalent linear model, 
doi: 10.1371/journal.pone.0080968.t004 



Table 5. Sensitivity of community indices to human 
pressures as defined and measured in EUROWATER, a 
special variant of CLC for freshwater common to the 
European scale. 



Index 


Value 


Meadow Farming Mix 


Urban 


Intens. Urban 


CSI 


Coef 


NS 


NS 


0.09 


0.04 


0.09 


R 2 = 0.27 


t-value 


NS 


NS 


10 


5 


9 




p-value 


NS 


NS 


<0.001 


O.001 


O.001 


LHT-COI 


Coef 


-0.008 


-0.01 


0.02 


NS 


0.01 


R 2 = 0.27 


t-value 


-3 


-4 


4 


NS 


3 




p-value 


0.003 


0.002 


<0.001 


NS 


0.01 


Niche-COl 


Coef 


-0.01 


-0.02 


-0.01 


-0.007 


-0.01 


R 2 = 0.27 


t-value 


-5 


-7 


-4 


-2 


-3 




p-value 


<0.001 


<0.001 


0.003 


0.03 


0.008 


Diet-COl 


Coef 


-0.01 


-0.009 


-0.01 


-0.02 


-0.02 


R 2 = 0.06 


t-value 


-7 


-4 


-4 


-9 


-6 




p-value 


<0.001 


<0.001 


<0.001 


<0.001 


<0.001 


Habitat-COl 


Coef 


NS 


-0.01 


0.05 


0.03 


0.05 


R 2 = 0.39 


t-value 


NS 


-3 


8 


5 


7 




p-value 


NS 


0.02 


<0.001 


O.001 


O.001 


Phylogenetic-COl 


Coef 


-0.03 


-0.03 


NS 


-0.02 


-0.03 


R 2 = 0.20 


t-value 


-4 


-4 


NS 


-3 


-3 




p-value 


<0.001 


<0.001 


NS 


0.009 


0.02 



(See table 4 for details.) 

doi: 10.1371/journal.pone.0080968.t005 



Here, we worked only with traits having a demonstrated 
functional role in fish biology. For example, morphological traits 
such as the mouth position or the length of barbell are linked to 
the diet and food acquisition [29,35]. Because our originality 
index is based on distance metrics we verified that it was not 
species richness dependant of the initial dataset. The life 
history traits index and the habitat-linked morphological index 
are less sensitive (R 2 >80, see S1) than the diet-linked 
morphological index and the niche index (respectively R 2 = 65 
and 68, see S1). Moreover, some community indices could be 
especially sensitive to one or a few species with extreme 
values of originality. In this study, it was the case with the 
European eel, which disproportionately increased the original 
community indices based on phylogeny and on habitat-linked 
morphology. This is the main purpose of this index, to weight 
unique species. However, because the European eel is 
classified as critically endangered at the national and global 
scales by the IUCN, these original indices also met in this 
particular case, the red list species indicator. Moreover, the 
European eel is a patrimonial species, there is a strong cultural 
heritage in France associated with this species for their fishing 
and cooking and thus for their taste, but also for their unique 
form and shape. For this last one based on the human vision, 
an originality index based on morphology could be common 
avenues for all species to be "objectively" quantified on the 
arguments develop by naturalists trying to preserve the unique 
forms and shapes that have been emerging on Earth. 
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Table 6. Sensitivity of community indices to human pressures as defined and measured in ONEMA, a special variant of CLC 
and EUROWATER for freshwater common to the French national scale. 



Index 


Value 


Meadow 


Farming 


Intens. Farming 


Mix 


Urban 


Intens. Urban 


CSI 


Coef 


-0.03 


NS 


NS 


0.08 


0.02 


0.07 


R 2 = 0.27 


t-value 


-5 


NS 


NS 


9 


3 


8 




p-value 


<0.001 


NS 


NS 


<0.001 


0.004 


<0.001 


LHT-COI 


Coef 


-0.03 


-0.02 


-0.02 


0.01 


NS 


NS 


R 2 = 0.28 


t-value 


-10 


-6 


-6 


3 


NS 


NS 




p-value 


<0.001 


<0.001 


<0.001 


0.02 


NS 


NS 


Niche-COl 


Coef 


-0.01 


-0.02 


-0.02 


-0.01 


NS 


-0.01 


R 2 = 0.28 


t-value 


-4 


-9 


-9 


-4 


NS 


-3 




p-value 


0.002 


<0.001 


<0.001 


0.001 


NS 


0.01 


Diet-COl 


Coef 


-0.008 


-0.01 


-0.007 


-0.008 


-0.02 


-0.01 


R 2 = 0.06 


t-value 


-4 


-6 


-3 


-3 


-8 


-5 




p-value 


<0.001 


<0.001 


<0.001 


0.006 


<0.001 


O.001 


Habitat-COl 


Coef 


NS 


NS 


-0.01 


0.05 


0.03 


0.05 


R 2 = 0.39 


t-value 


NS 


NS 


-3 


8 


6 


8 




p-value 


NS 


NS 


0.02 


<0.001 


<0.001 


<0.001 


Phylogenetic-COl 


Coef 


NS 


-0.04 


-0.03 


NS 


-0.02 


-0.02 


R 2 = 0.20 


t-value 


NS 


-6 


-4 


NS 


-2 


-2 




p-value 


NS 


<0.001 


<0.001 


NS 


NS 


0.02 



(See table 4 for details.) 

doi: 10.1371/joumal.pone.0080968.t006 



II. Spatio-temporal patterns 

We found that for the river fish communities, mapping each 
diversity component separately reveals partially congruent 
patterns between functional or phylogenetic originality, and 
specialization. All indices highlight Corsica, and to a lesser 
extent the Channel watersheds, as hotspots of originality and 
specialization. In contrast, the Seine watershed presents a less 
original and specialized fish community. For all the other 
watersheds, the different functional and phylogenetic COI and 
CSI are not congruent. For example, diet-linked and habitat- 
linked morphological COI present completely different patterns. 
The former highlights a strong originality in all small rivers 
especially mountainous streams, whereas the latter did not 
present any strong variation pattern at all, except for Corsica 
Island. 

These common patterns suggest that species occurring 
locally may be derived from regional species pools with similar 
biogeographical and evolutionary histories [63,64]. Moreover, 
for a given regional pool, species may respond to 
environmental gradients in different ways affecting the spatial 
distribution of the different biodiversity components and 
generating a spatial mismatch between functional and 
phylogenetic COI and CSI [64,65]. These results challenge the 
use of a single component as a surrogate for the others, and 
stress the need to first understand the different processes 
underlying each index and second to adopt a more integrative 



approach for conservation. One option to deal with the different 
messages given by the functional properties of communities 
and the resulting set of measurements is to be able to provide 
a hierarchy of their meanings depending on the context and 
perspectives or more reasonable to only use common patterns. 

A temporal effect has been detected for three original 
community indices (Table 3). Both morphological and Niche 
community indices significantly increased over time. This 
temporal dynamic should result from a global increase of the 
total species abundance, which has been shown (t= 5,09; 
p<0.001) and thus, may be the result of a global improvement 
of the entire river ecosystems. Indeed, global water quality has 
improved compared to the last century thanks to significant 
efforts to decrease organic substances [26,66]. Moreover, fish 
populations in Europe are still in their re-colonizing process 
since the last glaciations and some species are expected to 
extend their geographical area [67]. More precisely, each index 
is more influenced by the population dynamic of a few species 
presenting a high original value. The European eel has a very 
high original value especially for both habitat-linked 
morphological and niche indices, even though this species is 
declining (t= -6,53; p<0.001) which implies an important 
increase by other species in compensation. For the niche-COl, 
the global increase could be mainly linked with the population 
expansion of the European perch, northern pike, and minnow. 
The northern pike is very popular with anglers favoring their 
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Figure 6. Community indices response to land use artificialization with the CORINE Land Cover dataset. Colors from green 
to blue represent increasing pressures. 

doi: 10.1371/joumal.pone.0080968.g006 
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Figure 7. Niche-COl response to land use artificialization across watersheds. The Mix class is not represented because of 
limited sample size. 

doi: 10.1371/journal.pone.0080968.g007 
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introduction and thus, may have a positive impact on the 
population dynamic [26]. In the case of the diet-linked 
morphological-COl, it could be explained by three increasing 
species: the three spines stickleback, the common nase, and 
again the minnow. And in the case of the habitat-linked 
morphological-COl, the temporal increase could be mainly 
linked with the increase of two introduced species: the crucian 
carp and the pumpkinseed sunfish. Introduced species may 
increase faster than native species due to their rapid spread 
and their repeated introductions (accidental or deliberate). 

III. Community Indicators 

If originality and specialization approaches are completely 
different in their mathematical calculations, they describe 
complementary components of the functional properties of 
communities with similar expectations for their roles in 
ecosystem services. It has been theoretically and 
experimentally shown that the alteration of biodiversity disrupts 
ecological functions performed by species assemblages [68], 
and we know that species niche partitioning is fundamental in 
ecosystems properties [69]. Thus the more we are losing 
specialist or original species, the more we are losing 
irreplaceable functions in the ecosystem [6,70]. 

The theoretical background underlying the link between COI, 
CSI, and ecosystem functioning is growing in the literature, but 
it does not yet mean that both types of indices are relevant as 
functional biodiversity indicators. Community indicators have to 
be sensitive to anthropogenic pressures and give a clear and 
simple message to be technically and practically used by the 
targeted audience. Thus, the interpretation of indicators has to 
be as simple as possible and some communication qualities 
have to be accounting for. An indicator that is well recognized 
by all biodiversity stakeholders is more likely to be used in the 
future. According to the Millennium Ecosystem Assessment 
[17], one of the most important direct drivers of biodiversity loss 
and ecosystem service changes are land use modifications, 
including the physical modification of rivers. In this study we 
build innovative indices for community river fish in order to 
open the way for a new generation of indicators based on traits 
or niche linked with ecosystem functions. We compared each 
of our created community indices with the land use dataset 
(Figure 6). Interestingly the CSI present a very high score in 
urban areas. Fish communities are composed of specialist 
species in this artificialized habitat, probably because the 
environmental filter is very strong and species need to be 
specialists of this disturbed area such as the Black bullhead 
(Ameiurus melas), an invasive species which have a strong 
SSI. This pattern of urban specialist species in fish seems to be 
similar to the bird one, with urban specialist species such as 
pigeons (Colombia liva) or house sparrows (Passer 
domesticus) [21]. 

One COI appears particularly relevant to become a 
functional indicator of river fish communities, the COI based on 
species niche. This functional COI is sensitive to the different 
kinds of human disturbances with a simple interpretation: 
higher disturbance correlates with lower indicator values. 
Indeed, when we sorted land use on an artificialization gradient 
from natural to farming to urban areas, we observed a 



decrease in the niche-COl (Figure 6). Moreover, the three 
different variations of the CORINE Land cover data set give 
exactly the same results for the niche-COl, and thus we can be 
confident in its reproducibility over classification criteria. Finally, 
we evaluated its sensitivity to spatial scale. Is the pattern 
observed at the national scale still present at the regional 
scale? We observed a consistent pattern over the different 
watersheds, with the exception of the Mediterranean region 
(Figure 7), which may be due to the sample size. More 
interestingly, but not useful as an indicator, it may due to a 
local scale dependence or an eco-regional dependence. 
Indeed, Mediterranean watershed is very small and without any 
big river. In addition, all the urban areas are concentrated along 
the coast in this part of the Mediterranean eco-region. Further 
investigations need to be done to confirm the context of the use 
of this potential bio-indicator. 

Otherwise, Niche-COl as a functional biodiversity indicator 
encompasses species indicators like the CBD headline 
indicator "trend of selected species" because they consider a 
complete ecological group with their functions and their 
common dynamics. As a result, they carry more significant 
ecological information so that expectations and objectives for 
biodiversity stakeholders can be derived. Originality indices 
have already been used as indicators in biodiversity 
conservation contexts. Isaac et al. [71] built one called EDGE 
(Evolutionary Distinct and Globally Endangered) based on the 
phylogenetic originality and conservation status of species. At 
the community level, Mouillot et al. [70] have suggested 
functional and phylogenetic COI to evaluate conservation 
action areas. Because interspecific competition is more intense 
among species sharing common traits due to the limiting 
similarity principle [72], they expected that under protected 
areas competition might drive the better colonization or 
subsistence of the most original species because of niche 
complementary. We believe that even if one of the goals of this 
study was to develop functional biodiversity indicators for 
environmental policy makers, the Niche-COl could also be 
used at the scale of conservation reserves and may be used by 
managers of protected areas. 

Finally, we have to keep in mind that biodiversity indicators 
help to prioritize conservation actions to conserve ecological 
functions and in fine ecosystem "services". However, 
evaluation and measure alone are not sufficient in order to stop 
biodiversity loss, human pressures also must be limited. 
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